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5 ABSTRACT OF THE DISCLOSURE 

ABSTRACT 

A method for designing Wavelets for communications and 
radar which combines requirements for Wavelets and finite impulse 
response FIR filters including no excess bandwidth, linear 

10 performance metrics for passband, stopband, quadrature mirror 
filter QMF properties, intersymbol interference, and adjacent 
channel interference, polystatic filter design requirements, and 
non-linear metrics for bandwidth efficient modulation BEM and 
synthetic aperture radar SAR. Demonstrated linear design 

15 methodology finds the best design coordinates to minimize the 
weighted sum of the contributing LS error metrics for the 
respective performance requirements. Design coordinates are 
mapped into the optimum FIR symbol time response. Harmonic 
design coordinates provide multi-resolution properties and enable 

20 a single design to generate Wavelets for arbitrary parameters 
which include dilation, down-sampling, up-sampling, time 
translation, frequency translation, sample, rate, symbol rate, 
symbol length, and set of design harmonics. Non-linear 
applications introduce additional constraints. Performance 

25 examples are linear communications, BEM, and SAR. 

The present invention dcocribco a new method for the design 
e£ — multi - resolution — waveforms — fee — improve — system — performance, 

design — flexibility, applicability, a**d — reduce — the — costs — e# 

30 implementation, compared — fee — feke — current — a^4^ £ke — multi - 
resolution waveform designs arc a generalization of Wavelets to 
fe£e — Fourier domain. — This — generalization uses — design algorithms 
which include application metrics to improve system performance, 
allows a single design to be used for all scales, af*d allows 



1 



each dcoign to be characterized by a relatively few parameters. 

Preliminary simulations indicate that wavcf orm aftd filter 

performance — eaft — he — oignif icantly — improved — compared — te — current 
techniques — — communications — applications — fee — code — division 

5 multiple — access (CDMA) — wave forms — a**d — filters, — communications 

constant amplitude — bandwidth — efficient — waveforms, aftd — radar 

system applications. 

10 
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CLAIMS 

10 WHAT IS CLAIMED IS: 

Claim 1. (currently amended) An iterative eigenvalue least- 
squares LS method for designing digital mother Wavelets at 
baseband for A — means — fe* — the — design — — new — multi-resolution 
15 waveforms and filters , said method comprising steps: «i — the 
Fourier domain with properties which 
power spectral density PSD representative requirements for said 
mother Wavelet i|/ frequency o response V|/a(a)) in a multi- 
channel filter bank, specify 
20 a) passband frequency range for waveform transmission, 

b) stopband spacing between adjacent filters, 

c) bounds on ripple over said passband, 

d) stopband filter attenuation, 

e) rolloff with frequency outside stopband, 

25 f) quadrature mirror filters QMF require the sum of said 
PSD's for contiguous filter responses to be flat over 
deadband which is said stopband, 

g) symbol-to-symbol interference ISI, 

h) adjacent channel interference ACI, 

30 provide — extensions — e£ — the — Wavelet — concept — fee — the — Fourier 
domain or oquivalcntly the frequency domain 

provide single waveform — designs for all — of the waveform s at 

multiple scales 

35 
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provide — dcoigno — which — rcquirco — a — relatively — &ew — dcoign 
coordinates — compared to the — FIR time — rcoponoc samplco — £e*? 
the waveform 

provide dcoign methodologies which — eetn incorporate 

5 application metrics to improve the waveform performance 

provide — dcoigno — that — allow — feke — — e-£ — direct — dcoign 
mcthodologico — that — circumvent — fcke — need — te — solve — a Wavelet 
iterated filter bank construction to obtain the waveforms 

thereby providing improved flexibility fce meet fcke 

10 application goalo 

said LS error metrics to measure said requirements (a) -(h) are 

derived as functions of said Wavelet \(/(n) assuming 

i) T and 1/T are sample interval and sample rate equal to 

Nyquist sample rate, 

15 j) v|/ is real and symmetric about n=Q, 

k) n= 0, + /-1, . . . , +/-ML/2 digital index over said \y, 

1) M is interval between contiguous said vj/, 

m) 1/MT is said \\f symbol rate and channel-to-channel 

separation/ 

20 n) L is length of said in units of said M, 

said multiple-resolution properties require said LS metrics 

to be constructed as functions of said Wavelet Fourier 

harmonics iy k (k) with k=0, +/-!,. . + /-(N k -l) and 

N k >L is a design parameter, 

25 it is sufficient to use positive n=0,l,. . ML/2 and k=Q,l,. . 

. , N k -1 since said \\f (n) and nMk) are real and symmetric, 

ML/2+1 x N k matrix bw wherein "x" reads "by" maps \|/ k (k) into 

y/(n) to within a scale factor by equations 

Vj/(n) = E k bw(n+l / k+l)Vj/ k (k) for n>0, k>0, 

30 bw(n+l,k+l) = 1 for n=0, 

= 2 cos (27t nk/ML) otherwise, 

= row n+1, column k+1 element of bw, 

said M/(n) and \j/ k (k) are converted to column vectors by equations 

2 



ML/2 + 1 x 1 column vector of v|/(n) with elements 
h(l) = \|/(n=0) , 



h(n+l) = 2y(n) for n=l;2 y , . -ML/2, 

_h k = N k x 1 column vector of i|/ k (k) with elements 
h k (l) = V k (k=0) , 

h k (k+l) = 2y k (k) for k=l,2,. . . ,N k -l, 



said bw maps h k into h to within a scale factor by matrix equation 
h = bw h k/ 

LS error metrics for band = passband, stopband, and QMF deadband 

10 requirements are derived as quadratic forms in said h, 

said LS error metrics are converted by said bw mapping into 

quadratic forms in said h k equal to J(band)=h k 'R h k wherein 
said h k ' is the transpose of h k and said R is a real square 
symmetric matrix of LS errors in meeting said requirements, 
15 LS ISI,ACI error metrics J (ISI) , J (AC I) are derived as non-linear 
quadratic forms in h and converted by said bw matrix to the 
non-linear quadratic form in h k equal to J(ISI) =5e / 8e, 
J(ACI)=25E' 5E wherein 8E =AHh k is a column vector and 
matrix "A" in the matrix product AH is a function of said h 
20 hereby introducing said non-linearity, and said AH differ 

for ISI and ACI error metrics, 
LS cost function J is the weighted sum of said LS error metrics 

J = Z w(LS metric) J(LS metric) 

with summation over said LS metrics= passband, stopband, 
25 QMF deadband, ISI, ACI with normalized weights 

Zw(LS metric)-!, 
said weights, are free design parameters, 

said iterative eigenvalue LS algorithm at each step finds the 

optimum eigenvalue and eigenvector which minimize said 

30 quadratic form J in h k for a constant said "A", 

said eigenvector is the optimum h k which minimizes said J and 

said bw equation derives the corresponding optimum h which 

minimizes said J, 
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step 1 in said iterative algorithm finds said optimum eigenvalue/ 
eigenvector, h k , h of J reduced by deleting said non-linear 
ISI and ACI LS quadratic error metrics, 
said h is used to evaluate said "A" matrices for step 2, 
5 step 2 finds said optimum eigenvalue, eigenvector, h k , h for 

minimum J using said "A" from step 1, 

said h is used to evaluate said "A" for step 3, 

steps 3,4, etc. continue until said minimum J converges to a 

steady value and, 

10 said optimum \|/ k (k) uses said bw to calculate optimum v|/(n) for 

implementation as said Wavelet FIR digital waveform and 

filter time response. 



15 

Claim 2. (currently amended) An LS method for designing digital 
mother Wavelets at baseband for multi-resolution waveforms and 
filters, 2. A means for tho dcoign of multi rcoolution wavoformo 
and filters in the — frequency domain or cquivalcntly tho — Fourier 
20 e£ — fcke — fc— £ — space — fef — waveform — af*3 — filter — applications, with 
propcrtico which said method comprising steps: 

said PSD waveform representative requirements and assumptions 
are recited in (a)-(n) in Claim 1, 

said multiple-resolution properties require said LS metrics 

25 to be constructed as functions of said iy k (k), 

said LS error metrics for said passband, stopband, and QMF 

deadband requirements are derived as squared vector norm 

functions of said h and converted by said bw matrix into 



J (band) = Bh 
II 



- wherein 



Bh is the vector norm of the 
k|| 



30 column vector Bh k and said B is the matrix of LS errors in 

meeting said requirements and wherein said squared vector 

norm is suitable for LS optimization, 
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provide a means for single dcoign to be used for all of the 
waveforms at — the multiple resolutions or frequency bands 

provide — a moans — for designs — that — require — a — relatively few 
design — coordinates — compared — fee — the — number — en§ — digital 
5 samples covered by the Wavelet 

provide a means for the designs to include frequency 

and — time — application — metrics — fee — improve — fehe — waveform 
performance 

LS ISI,ACI error metrics J ( ISI ) , J (ACI ) are derived as squared 

1 0 vector norm functions equal to J(ISI = jjs e|| 2 , J (ACI) =2 1|6 e|| 2 

using said column vectors 5E =AHh k in claim 1, 
LS cost function J is said weighted sum of said LS error metrics 

equal to J=£w(LS metric) J (LS metric) defined in claim 1, 

an LS gradient search algorithm finds optimum h k (k) to 

15 minimize J, 

step 1 of said LS gradient search algorithm uses a Remez- 

exchange algorithm to find said optimum h k (k) for said J 

reduced to said passband and stopband LS metrics, 

20 provide a means — to use direct design methodologies — in the 

frequency time — domain — e — fe^e — arti — e# — feke — waveforms — and 
filters for the multiple frequency bands 
step 2 uses the estimated h k (k) from step 1 to initialize said 

gradient search, 

25 step 3 selects one of several available gradient search 

algorithms, gradient search parameters, and stopping rules, 

step 4 implements said algorithm, parameters, and stopping rule 
selected in step 3 to derive said optimum h k (k) to minimize 
J and, 

30 said optimum h k uses said bw to calculate optimum v|/(n) for 

implementation as said Wavelet FIR digital waveform and 
filter time response. 
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Claim 3. (currently amended) Wherein said mother Wavelet in 
Claims 1,2 generates multi-resolution dilated Wavelets, 
comprising steps and design: 
5 said Wavelet parameters are 

a) scaling parameter p dilates sampling by factor 2 A p 
equivalent to sub-sampling by factor 2 A p, 

b) translation parameter q translates said y; by qM digital 
samples, 

10 c) frequency offset k is set by design, 

d) symbol repetition interval said M remains constant, 

e) Wavelet length said L in units of said M remains 
constant, 

step 1 uses said design harmonics \y k to generate said FIR time 

15 response \|/(n p ) at baseband with said bw equation 

V(n p ) = Z k bw(n p +l, k+1) \\f k (k) 

recited in claim 1 with 



bw = 2 cos(2n n p k/ML) for n p >0 
wherein n p =n/2 A p is n sub-sampled or equivalently dilated 

20 by the factor 2 A p, 

step 2 uses said y/(n p ) to construct said multi-resolution 
Wavelet i|/ p , q ,M,L,k with equation 

^,q,M,L,k = 2 A (-p/2) \i;(n p -qM) exp(j27t kn p /ML) 

which is said FIR time response for parameters p f q f M f L f k 

25 wherein the subset p,M,L are the scale parameters, 

ft — method — fe^ — the — design — e£ — multi - resolution — waveforms — which 

allows — Fourier — domain — techniques — fee — be — used. Properties — eeft 

include some or more of the listed properties in — (^H — aftd — (-2-)- 
30 design of said multi-resolution Wavelet includes 

f) said T for n is increased to T2 A p for n p , 

g) said 1/T is reduced to l/T2 A p, 
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h) said \\f symbol rate 1/MT equal to said channel-to-channel 
separation is reduced to l/MT2 A p in Hz and, 

i) said \\f length (ML+1)T in seconds is stretched to 
(ML+l)T2 A p in seconds. 



Claim 4. (currently amended) Wherein said mother Wavelet in 
claims 1,2 generates multi-resolution constant sample rate 
10 dilated Wavelets, comprising steps and design: 
said Wavelet parameters are 

a) said p dilates said \\f to increase said length from 
ML+1 to MpL+1 where M p =M2 A p is the dilated interval 
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b) 


said 


q 


translates said vj/ by qM p digital samples, 


c) 


said 


k 


is set by design, 


d) 


said 


M p 


=M2 A p is dilated M, 


e) 


said 


L 


remains constant , 



step 1 uses said design harmonics V|/ k to generate said FIR time 

20 response \y(n p ) at baseband with said bw equation 

V(n) = Zk bw(n+l,k+l) v|/ k (k) 

recited in claim 1 with 



bw = 2 cos (271 nk/M p L) for n>0, 



step 2 uses said \\f (n) to construct said multi-resolution Wavelet 

25 ^p, q ,M,L,k with equation 

M>P,q,M,L,k = 2 A (-p/2) \j/(n-qM p ) exp(j27i kn/M p L) 

which is said FIR time response for parameters p,q,M,L,k, 

design of said multi-resolution Wavelet includes 
f) said T remains constant, 

30 g) said 1/T remains constant, 

h) said \y symbol rate 1/MT equal to said channel-to-channel 
separation is reduced to 1/M P T=1/MT2 A p in Hz and, 

i) said \y length (ML+1)T in seconds is stretched to 
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(ML2 A p+l)T in seconds. 
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A method for the dooign of multi rcoolution waveformo which 
eaR — incorporate — Fourier — domain — tcchniquoo — into — dcoign 
methodologies — which — eae — include — analytical — aftd — iterated 
filter bank construction design tcchniquoo. 

10 

Claim 5. (currently amended) Wherein said mother Wavelet in 
claims 1,2 generates multi-resolution up-sampled Wavelets, 
15 comprising steps and design: 
said Wavelet parameters are 

a) said p up-samples said digital sampling rate 1/T to 2 A p/T, 

b) said q translates said \\f by qM digital samples, 
c) said k is a design parameter, 

20 d) said M is constant 
e) said L is constant 
step 1 uses said design harmonics h fL vi/ k to generate said FIR 

time response v|/(n p ) at baseband with said equation 

V(n p ) = X k bw(n p +l,,k+l) v k (k) 

25 recited in claim 1 with 

bw = 2 cos (271 n p k/ML) for n p >0 
wherein n p is n up-sampled by the factor 2 A p and defined 
by equations 

n p = n p + n 2 A p 
30 n p = 0,1,2, . . .,2 A p-l 

wherein n p is the index over the additional samples added 

to each sample n by said up-sampling, 

step 2 uses said ij/(n p ) to construct said multi-resolution Wavelet 
vj£p, q ,M,L,k with equation 
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V£>,q,M,L,k = 2 A (-p/2) y(n p -qM) exp(j27i kn p /ML) 

which is said FIR time response for parameters p,q,M,L,k, 



A method for the analysis and design of multi rcoolution 
waveforms using Fourier domain techniques which take 
advantage of the new invention disclosures on the 
characterization and design of multi resolution waveforms 
10 in the — Fourier domain. 

design of said multi-resolution Wavelet includes 

f) said T is decreased to T/2 A p y 

g) said 1/T is increased to 2 A p/T, 

h) said \y symbol rate 1/MT equal to said channel-to-channel 
15 separation is increased to 2 A p/MT in Hz and, 

i) said \\f length (ML+1)T in seconds is reduced to 
(ML +l)T/2 A p in seconds. 
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Claim 6. (currently amended) Wherein said multi-resolution 
Wavelets in Claims 1-5 have properties comprising: 
said scale parameters p,M,L and said design parameter 1/T specify 

25 said multi- resolution Wavelets at baseband and said q, k 

specif ify time, frequency translations from baseband, 

said design harmonics i|/k(k) of mother Wavelet are said design 

coordinates for multi-resolution Wavelets, 

said design harmonics yMk) use said bw matrix to generate said 
30 multi-resolution Wavelet baseband time response \y(n) for 

dilation, dilation of Wavelet length, and up-sampling as 
recited in Claims 3-5 and which is translated in time and 
frequency to said multi-resolution Wavelet i|/ p , q , M , l, k±_ 



said design harmonics i|/ k (k) are few in number compared to said 
V (n) , 

said \\f is designed to support a bandwidth-time product B f T=l+a 

with no zero excess bandwidth g=0, 

5 said multi-resolution Wavelets are designed to behave like an 

accordion in that at different scales said Wavelets are 

stretched and compressed versions of the mother Wavelet with 

appropriate time and frequency translation, 

10 A — aew — formulation — i-fi — f5-) — fe* — multi - resolution — wavof orm — as — a 
function — — feke — de — mult i - rcoolut ion — waveform — which — adds — feke 
conccpt — e£ — a — frequency — index — that — allows — the — multi resolution 
waveform — to — be — placed — arbitrarily — throughout — fete — fe— £ — space 
thereby — — avoiding — fefee — restrictions — e£ — the — Wavelet — iterated 
15 filter construction for tiling a t - f space, — and 2) — allowing the 
new multi - rcsoloution waveforms — to be — used — for multi - resolution 

communications aed fe^ — bandwidth on - demand — communications 

application, — in place of traditional Wavelets. 
said optimization techniques in claims 1,2 assume said \|/(n) 
20 symmetric about n=0 and are applicable to other arrangements 

of \|/(n) with self-evident modifications, 
optimization algorithms for finding said optimum set of \\f k (k) 

use linear LS waveform and filter design methods recited in 
claims 1,2 and also use other methods and, 
25 said linear waveform and filter LS design methods can be modified 

to design waveforms for applications including bandwidth 

efficient modulation BEM and synthetic aperture radar SAR. 
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BACKGROUND OF THE INVENTION 



I, Field of the Invention 

10 TECHNICAL FIELD 

The present invention relates to CDMA (Code Division 
Multiple Access) cellular telephone and wireless data 
communications with data rates up to multiple Tl (1.544 Mbps) and 
higher (>100 Mbps), and to optical CDMA. Applications are 
15 mobile, point-to-point and satellite communication networks. More 
specifically the present invention relates to a new and novel 
means — for a new approach to the design of waveforms and filters 
using mathematical formulations which generalize the Wavelet 
concept to communications and radar. 
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BACKGROUND ART 



II. Description of the Related Art 



Multi-resolution waveforms used for signaling and/or 
5 filters and which are addressed in this invention are defined to 
be waveforms of finite extent in time and frequency, with scale 
and shift properties of multi-resolution over the time-frequency 
(t-f) space. These waveforms can also be referred to as multi- 
scale waveforms, and include multi-rate filtering and the 
10 Wavelet as special cases. The emphasis will be on digital 
design and applications with the understanding that these multi- 
resolution waveforms are equally applicable to analog design and 
applications . 

15 Background art consists of the collection of waveform and 

filtering design techniques which can be grouped into six broad 
categories. These categories are: CI) least squares (LS) design 
algorithms for filters and waveforms that design to 
specifications on their frequency response, C2) analytic 

20 filters and waveforms which are specified by a few free design 
parameters that can be sub-categorized into current applications 
and primarily theorectical studies, C3) combinations of CI and 
C2 for greater flexibility in meeting communications and radar 
performance goals, C4) special design techniques to control the 

25 noise levels from intersymbol interference (ISI) and adjacent 
channel interference (ACI) in the presence of timing offsets for 
multiple channel applications, C5) Wavelet filter design using 
scaling functions (iterated filter banks) as the set of design 
coordinates or basis functions, C6) filter and waveform design 

30 techniques for non-linear channels and in particular for 
operation in the non-linear and saturation regions of a high 
power amplifier (HPA) such as a traveling wave tube (TWT) or a 
solid state amplifier, and C7) LS dynamic filters derived from 
discrete filtering and tracking algorithms that include adaptive 

35 equalization for communications, adaptive antenna filters, 
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Wiener filters, Kalman filters, and stochastic optimization 
filters . 

Category CI common examples of LS digital filter design are 
5 the eigenvalue algorithm in "A New Approach to Least-Squares FIR 
Filter Design and Applications Including Nyquist Filters" [7] and 
the Remez-Exchange algorithm in "A Computer Program for Designing 

Optimum FIR Linear Phase Filters" . -f8-b The eigenvalue algorithm 

is a direct LS minimization and the Remez-Exchange can be 

10 reformulated as an equivalent LS gradient problem through proper 
choice of the cost function. Both LS algorithms use the FIR 
(finite impulse response) digital samples as the set of design 
coordinates. LS design metrics are the error residuals in 

meeting their passband and stopband ideal performance as shown 

15 in FIG. 3. Category C2 common examples of analytical waveforms 
and filters for system applications are the analog Chebyshev, 
Elliptic, Butterworth, and the digital raised-cosine, and square- 
root raised-cosine. For theorectical studies, common examples 
are polyphase multirate filters, quadrature mirror filters 

20 (QMF) , and perfect reconstruction filters. Although these 
theorectical studies have yet to yield realizable useful filters 
for system applications, their importance for this invention lies 
in their identification and application of ideal performance 
metrics for filter designs. Category C3 common digital example 

25 is to start with the derivation of a Remez-Exchange finite 
impulse response — (-FIR-}- filter and then up-sample and filter with 
another bandwidth limiting filter. This results in an FIR over 
the desired frequency band that is larger than available with the 
Remez-Exchange algorithm and with sidelobes that now drop off 

30 with frequency compared to the flat sidelobes of the original 
Remez-Exchange FIR. A category C4 common example is to select 
the free parameters of the category C3 filter in order to 
minimize the signal to noise power ratio of the data symbol (SNR) 
losses from the ISI, ACI, and the non-ideal demodulation. A 

35 second common example is start with a truncated pulse whose 
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length is short enough compared to the symbol repetition 
interval , to accommodate the timing offsets without significant 
impact on the ISI and ACI SNR losses. This shortened pulse can 
then be shaped in the frequency domain. 

5 

Category C5 Wavelet filter design techniques discussed in 
the next section will serve as a useful reference in the 
disclosure of this invention. Category C6 common example is the 
Gaussian minimum shift keying (GMSK) waveform. This is a 

10 constant amplitude phase encoded Gaussian waveform which has no 
sidelobe re-growth through a non-linear or saturating high power 
amplifier JHPA)_. Category C7 common examples are the adaptive 
equalization filter for communication channels, the adaptive 
antenna filter, and the Kalman filter for applications including 

15 target tracking and prediction as well as for equalization and 
adaptive antennas . 

Minimizing excess bandwith in the waveform and filter 
design is a key goal in the application of the C1,...,C6 design 

20 techniques for communications and radar. Excess bandwith is 
identified as the symbol a in the bandwidth-time product 
BT s =l+a where the two-sided available frequency band is B and the 
symbol repetition interval is T s . Current performance capability 
is represented by the use of the square-root raised cosine (sq-rt 

25 rc) waveform with ct= 0.22 to 0.4 as shown in FIG. 46. The goal 
is to design a waveform with a=0 within the performance 
constraints of ISI, ACI, passband, sideband, and passband 
ripple. This goal of eliminating the excess bandwidth corresponds 
to the symbol rate equal to the available frequency band 1/T S = 

30 B/(l+a) = B for a=0 . This symbol rate 1/T S =B is well known to be 
the maximum possible rate for which orthogonality between 
symbols is maintained. A fundamental performance characteristic 
of our new waveform designs is the ability to eliminate excess 
bandwidth for many applications. 
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Scope of this invention will include all of the waveform 
and filter categories with the exception of category C4 special 
design techniques and category C7 LS dynamic filters. Emphasis 
5 will be on the category C5 Wavelets to establish the background 
art since Wavelets are multi-resolution waveforms that can 
eliminate the excess bandwidth and have known design algorithms 
for FIR waveforms and filters. However, they do not have a design 
mechanism that allows direct control of the ISI, ACI, passband, 
10 sideband, and passband ripple. Category C2 theorectical studies 
also eliminate the excess bandwidth. However, they are not multi- 
resolution waveforms and do not have realizable FIR design 
algorithms. So the emphasis in background art will be on 
Wavelets whose relevant properties we briefly review. 



Wavelet background art relevant to this invention consists 
of the discrete Waveform equations and basic properties, 
application of Wavelets to cover a discrete digital time- 
frequency (t-f) signal space, and the design of Wavelets using 

20 the iterated filter construction. Wavelets are waveforms of 
finite extent in time (t) and frequency (f) over the t-f space, 
with multi-resolution, scaling, and translation properties. 
Wavelets over the analog and digital t-f spaces respectively are 
defined by equations (1) and (2) as per Daubechies' s "Ten 

25 Lectures on Wavelets", Philadelphia: SIAM, 1992 
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rof crcncG — [1] . 



Continuous wavelet 



/*\ I r 1/2 



a 



(1) 
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Discrete wavelet 



e . j j-i/2 f n-b 
a 



(2) 
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where the two index parameters M a,b" are the Wavelet dilation and 
translation respectively or equivalently are the scale and shift. 
The \|/ is the "mother" wavelet and is a real and symmetric 
localized function in the t-f space used to generate the doubly 

i r 1/2 

5 indexed Wavelet i|/ ab - The scale factor "|a| " has been chosen to 

keep the norm of the Wavelet invariant under the parameter change 
"a,b". Norm is the square root of the energy of the Wavelet 
response. The Wavelets \|/ a#b and y are localized functions in the 
t-f space which means that both their time and frequency lengths 
10 are bounded. The discrete Wavelet has the time "t" replaced by 
the equivalent digital sample number "n" assuming the waveform is 
uniformly sampled at W T" second intervals. 

Wavelets in digital t-f space have an orthogonal basis that 
15 is obtained by restricting the choice of the parameters "a,b" to 
the values a=2" p , b=qM2 p where A p,q" are the new scale and 
translation parameters and "M" is the spacing or repetition 
interval T S =MT of the Wavelets (which from a communications 
viewpoint are symbols) at the same scale "p" . Wavelets at "p,q" 
20 are related to the mother Wavelet by the equation +±+ 

^ pq (n) = 2-^K2- p n-qM) (3) 

25 

where the mother Wavelet is a real and even function of the 
sample coordinates. The orthonormality property means that these 
Wavelets satisfy the orthogonality equation with a correlation 
value equal to "1". 

30 
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Z^P-q^m = 1 iff both p = k and q = m 

n 

= 0 otherwise 



Wavelet representation of a digital t-f space starts with 
selecting an N sample time window of a uniform stream of 
digital samples at the rate of 1/T Hz (1/second) equivalent to a 
5 W T" second sampling interval. The N point or sample t-f space 
in FIG. 1 illustrates a Wavelet representation or "tiling" with 
Wavelets that are designed analytically or by an iterated filter 
construction. 

10 The t-f space in FIG. 1 is partitioned or covered or tiled 

by a set of Wavelet subspaces §V p ,p = 0 5 l 5 ...,/w-l} 1 where N=2 m . Each 
Wavelet subspace W p 1 at scale "p" consists of the set of Wavelet 
time translations {q = 0,1,..., N/2 P+1 -1} over this subspace. These 
Wavelet subspaces are mutually orthogonal and the Wavelets within 

15 each subspace are mutually orthogonal with respect to the time 
translates. This N-point t-f space extends over the time 
interval from 0 to (N-l)T 2 where — T — is — fefee — digital — oampling 
interval, — and over the frequency interval from 0 to (N-l) 3 in 
units of the normalized frequency f NT 4. 

20 

The iterated filter bank in FIG. 2 is used to generate the 
Wavelets which cover the t-f space in FIG. 1. Each filter stage 
5 consists of a high pass filter (HPF) and a low pass filter 
(LPF) . Output 6 of the LPF is subsampled by 2 which is 
25 equivalent to decimation by 2. This t-f space space is an N- 
dimensional complex vector metric space VI. At stage m in the 
iterated filter bank, the remaining t-f space V m _ 2 8 is 
partitioned into V m _ { 9 and the Wavelet subspace W m _i 10 . 

30 Scaling functions and Wavelets at each stage of this filter 

bank satisfy the following equations 
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Kn)=r ,/2 Sh q ^(2n-q) Vq 

(5) 



^(n)=2- ,/2 Sg q ^(2n-q) Vq 

4 



5 where cp is the scaling function, \\f is the Wavelet, HPF P 

coefficients are { h q9 Vq } , LPF P coefficients are {g 9 ,Vq}, and 

the equations apply to the stages 0, 1, ...,m-l . Identifying the 
scale parameter and using the previous Wavelet formulations 
enable these equations to be rewritten for stages p=0, l,...,m-l 
10 as 



y p = Zh q <p pA q Vp 

(6) 

<P P = Sg q <p p - hq Vp 

q 



For our application the HPF P and LPF P are quadrature mirror 
15 filters (QMF) with perfect reconstruction. This means they cover 
the subspace V p with flat responses over the subband frequency 
including the edges of the frequency subband, and the HPF P 
coefficients are the frequency translated coefficients for the 
LPF P : {g q =(-l) q h q , Vq}. 

20 

Wavelet design using iterated filter bank starts with the 
selection of the scaling functions. Starting with a primative 
scaling function such as the one proposed by Daubechies [1], — one 
can use the iterated filter construction given by equations (5) 
25 and (6) to derive successive approximations to a desired scaling 
function which has properties that have been designed into it by 
the selection of the filter coefficients {g^Vq} at each level of 
iteration. The Wavelets can be derived from these scaling 
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functions using the iterated filter construction or scaling 
equations (5) and (6) . 

Another use of the iterated filter construction is to 
5 design the scaling functions as Wavelets thereupon ending up with 
a larger set of Wavelets for multi-resolution analysis and 
synthesis as illustrated by Go if man ' s Wavelets in "Wavelet 
analysis and signal processing". — Design examples abound in 
the open literature and the recent publications — hH-y — \r^h — [5] , 
10 -£-€-} — arc representative illustrations of the diverse applications 
of Wavelets. 

SUMMARY OF INVENTION 

15 

SUMMARY OF THE INVENTION 

This invention is a new approach to for the design of new 

20 Wavelet m ulti -resolution waveforms that improves their 

performance for engineering and scientific applications, and in 
particular for applications to communications and radar. Current 
practice is to 1) design traditional multi-resolution waveforms 
in the time domain using metrics which specify frequency 

25 performance and without consideration of Wavelet propertiesT aftd 

or to 2) design Wavelet multi-resolution waveforms using the 
time-domain iterated multi-resolution filtering approach with a 
set of scaling functions used to perform the filtering and 
without direct considerations of the frequency performance. 
30 These two separate design approaches yield fundamentally 
different waveforms. Our ncw This invention provides a means to 
combine these two approaches to generate a new multi-resolution 
waveform with the best properties of the traditional multi- 
resolution waveforms and the Wavelet multi-resolution waveform. 
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Our Yhis invention introduces t&e — new innovations for the 
design of our new waveform that 1) include the ISI (inter symbol 
interference) , ACI (adjacent channel interference) , and QMF 
5 (quadrature mirror filter) requirements into the design algorithm 
along with non-linear modifications to the traditional passband 
and stopband frequency requirements which are used in the Remez- 
Exchange and eigenvalue LS algorithms — ifi — referenced — [-84-7 — 
rcopcctivcly , 2) select the frequency harmonics as the design 
10 coordinates, 3) use an appropriate subset of the available 
Fourier domain frequency harmonics as the design coordinates with 
the property that this subset of coordinates is a basis for the 
multi-resolution waveform design, and 4) design the multi- 
resolution waveform for no excess bandwidth a=0 . 

15 

These innovations support the development of our ftew 
waveforms 1) that are generalizations of Wavelets in the 
frequency domain using a means which makes them useful for 
communications and radar over their t-f space, 2) that can be 

20 designed to be an orthonormal basis or orthonormal set of 
coordinates over the available frequency bandwidth with the 
implicit property that the excess bandwidth vanishes a=0— to 
within the accuracies allowed by communication and radar design 
implementations, and 3) that have multi-scale properties which 

25 allow a single dc waveform design to be used to uniquely define 
the complete set of new waveforms over the t-f space for multi- 
scale applications^ Specific design algorithm examples developed 
in this invention disclosure are the LS design algorithms that 
use eigenvalue and gradient search techniques in the frequency 

30 domain to find the best waveform design which minimizes the 
corresponding LS error residual cost function that is a weighted 
linear sum of the residual error metrics. 
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Waveform performance calculated in this invension 
disclosure from the application of these LS design algorithms, 
demonstrates the capability of this invention to provide a means 
to design waveforms that are improved over current practice. 
5 There are other algorithmic design algorithms which can be 
realized by modifications to this invention. An example given 
in this invention disclosure is the modification of these LS 
algorithms for application to constant amplitude bandwidth 
efficient (BEM) waveforms for communications. 

10 

BRIEF DESCRIPTION OF THE DRAWINGS AND 
THE PERFORMANCE DATA 

BRIEF DESCRIPTION OF DRAWINGS AND PERFORMANCE DATA 

15 The above-mentioned and other features, objects , design 

algorithms, and performance advantages of the present invention 
will become more apparent from the detailed description set forth 
below when taken in conjunction with the drawings and performance 
data wherein like reference characters and numerals denote like 

20 elements, and in which: 

FIG. 1 is a Wavelet aa— N-point t-f space extending over the 
time interval (0, (N-l)T] and the frequency interval (0, (N- 
1)/NT], which is tiled or covered by a set of orthonormal 
25 Wavelets at the scales p=0, 1, ...,m-l . 

FIG. 2 is a Wavelet iterated filter bank used to generate 
the set of Wavelets which tile or cover the t-f space in FIG. 1. 

30 FIG. 3 illustrates the power spectral density (PSD) of our 

new waveform at dc, and the set of stopband and passband design 
requirements . 
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FIG. 4 is the flow diagram of the LS metrics and cost 
functions and the final cost function used to find the optimal LS 
solution for our waveform. 

5 FIG* 5 is flow diagram for the LS recursive solution 

algorithm to find the optimal harmonic coordinates that define 
the LS solution for our waveform. 

FIG. 6 plots the dc PSD in dB units for the — our new 

10 waveforms vs. the normalized frequency fT s for t-he our new 

waveform designed with the LS algorithms and for the sq-rt r-c 
waveforms with excess bandwidth parameter a= 0.22, 0.40. 

15 

FIG. 7 plots the dc PSD in dB units for our wavoformo vs. 
the normalized frequency fT s for the our new waveform designed 
with modifications to the LS algorithms for application to 
constant amplitude BEM communications, and for a Gaussian 
20 minimum shift keying (GMSK) waveform. 

FIG. 8 plots the amplitude of the dc radar ambiguity 
function vs. the normalized frequency fT p and normalized time tT c 
25 where the pulse time T p and chip time T c are both equal to the 
communications symbol interval T s for this example, for the our 
new waveform and an unweighted chirp waveform. 

30 



35 
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DISCLO S URE OF INVENTION 

DISCLOSURE OF THE INVENTION 



5 

New waveforms in this invention disclosure are 
generalizations of Wavelets in t-f space which enable them to be 
useful for communications and radar applications. This 

10 generalization is accomplished by 1) the addition introduction 
of a frequency translation, by 2) relaxing — changing the 
orthonormality condition in equation (4) to apply to waveforms 
within the same space {q} and over the scales {p} with the 
inclusion of the frequency translation, and 3) by their 

15 characterization and design in the Fourier domain. 

With frequency translation the analytical formulation of 
these new waveforms as a function of the baseband or mother 
waveform centered at dc_(— dc refers to the origin f=0 of the 
20 frequency space) becomes 



New waveform ao a function of dc waveform 

V M ,00 = 2- p V(2- p n -qM)e i2rf ' (p - r)nT (7) 

25 

where f c {p,r) is the center frequency of the frequency translated 
dc waveform, at scale "p" and frequency index V. The purpose 
of the frequency index "r" is to identify the center frequencies 
30 of the waveforms at the scale "p" in the t-f space. The dc or 
baocband or mother waveform io the generalization of the mother 
Wavelet for multi rcoolution waveforms. 
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These waveforms satisfy the complex orthonormality 
equations 



Orthogonality cquationo 

Y\if = 1 iff p = kandq - m and r = v 

£ y p,q,r Y k,m 9 v 

= 0 otherwise (8) 



where is conjugation, and are generalizations of the 

orthonormality equations for the analytical Wavelets in (4) . 

10 fee — new — waveforms — in — equation — W expand — the — Wavelet 

analytical formulation fee include a frequency — variable. 

Wavelets arc — functions of the ocalc and translation parameters 
"p,q". The concept of an additional parameter to provide an added 
degree of flexibility in their tabulation was first introduced by 

15 the — Coifman library of Wavelets — f-2-] — which use — an — oocillaton 

parameter tabulation that roughly corresponds fee t^e 

frequency of oscillation of the Wavelets. The frequency variable 

together with the Fourier domain design arc entirely new means 

£&e — deriving — these — new — waveforms — as — generalisation — e£ — the 

20 traditional — Wavclcto — and — their — modification — with — the — added 
parameter in — [2] . 

New — Our new waveforms are generalizations of Wavelets in 
the frequency domain using a means which makes them useful for 

25 communications and radar in the t-f space. The basis vectors for 
this metric space V consist of a subset of the admissible set 
of scaled and translated waveforms {^ Mr ,Vp,q,r} derived from the 
dc waveform y/ as per equation (7) . An admissible waveform is 
any combination that covers F=t-f space. We are interested in 

30 the Fourier domain representation of the dc waveform y/ in V , 
and in particular in a subset of the discrete Fourier transform 
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(DFT) harmonic coefficients over the Fourier domain which we 
intend to use as the design coordinates. Starting with the z- 
transform and continuous Fourier transform, the DFT harmonic 
coefficients are defined by the following equations 



DFT Harmonic Coefficients (9) 
^y(n) z* z-transform 

n 

£y(n) e _uan Fourier transform 

n 

^y(n) DFT harmonic coefficients for V£ 

n 

k n 

sin((<y/2- t&JNW)! Nsin(6V2 -;zk/N') 

k 

^y/ k [Harmonic interpolation for "k"] 

k 

_i2Tr2TTkn/ 
= e 

{k} = DFT frequency or harmonic coefficients such that f k =k/ISTT 

where f k is the harmonic frequency corresponding to "k" 
N 1 = length of qj(n) 

where ^{co) -\y/{co)f is the power spectral density of the dc 
waveform. These equations define the frequency representation of 
our new waveforms in terms of the available set of harmonic 
coefficients, which set is considerable larger than required for 
most applications. 

Mew — Our new waveforms are an orthonormal basis with no 

excess bandwidthT arc properties — that — we — will — demonstrate . . 

Orthomality and no excess bandwidth £H?e — properties which are 



y/{z) = 
y/(co) = 

Vk = 

where 
y/((D) 
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asymptoticly approached by our — new wave forms — to within design 
accuracies inherent in communications and radar. To proceed we 
need to identify the structure of the dc waveform in V . We 
start with definitions for the parameters and coordinates in the 
5 following equation (10) , The waveforms derived from the dc 
waveform will be designed to be orthogonal over both time 
translates "MT" and frequency translates XX 1/L£MT" which 
respectively correspond to the Wavelet symbol spacing T S =MT and 
the adjacent channel spacing 1/T S . This means the orthogonal 

10 spacing of the waveforms in V are at the time-frequency 
increments (MT, 1/i^MT) = (T s . 1/T S ) . In the interests of 
constructing our f*ew— orthonormal multi-resolution waveforms to 
cover V it will be convenient to" assume that M,L are powers of 
2. We need the following definitions for the parameters and 

15 coordinates. 



Parameters and Coordinates (10) 

20 
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N' = Length of \|/ which is an even function about the center and 
which spans an odd number of points or samples 
= ML + 1 where M, L are assumed to be even functions for 

convenience of this analysis 
= Number of points of y 
M - Sampling interval for y/ 

- Spacing of y/ for orthogonality 
L = Length of \|/ in units of the sample interval M 

= Stretching of xf over L sample intervals 
n = n 0 + n,M 

= partitioning into an index n 0 over the sample 

length n 0 = 0,1,...,M- land an index n, = 0,1...,L-1 
over the sample intervals 
k =k 0 + k,L 

= partitioning into an index k 0 over the harmonic 
frequencies k 0 = 0,1,...,L-1 corresponding to the 
stretching and an index k x = 0,1,...,M over the harmonics 
frequencies corresponding to the admissible frequency 
slots for if/ 

The harmonic design coordinates are selected using the 
following observation. For most applications and in the following 
development it is assumed that the waveforms are spectrally 
5 contained in the frequency interval 1/L3- MT corresponding to the 
frequency spacing. This suggests the harmonic design coordinates 
be restricted to the subset of L harmonics {k 0 =0 / 1, . .., L-l} 
covering this spacing. These L harmonics correspond to the 
stretching of the mother waveform over the L repetition 
10 intervals. 

Obviously, for some applications as will be demonstrated 
later, the spectral containment is spread out over several 1/LT 
frequency increments whereupon one must increase the subset of 
15 design harmonics to possibly 2L, 3L or larger. 
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The DFT equations for the dc waveform in (9) when rewritten 
in terms of the L harmonic design coordinates {^,Vk 0 } become: 

5 

DFT equations for dc waveform 

y/ k = ^y(n) W^° n harmonic design coordinates 

n 

y/(n) =(1/N') E*F k W*° n new waveform defined in terms (11) 
10 of the L harmonic design coordinates {y^Vko} 

the use of these L harmonic design coordinates is sufficient 
under time translates to be a basis for the corresponding 

15 subspace of V . This means these waveforms provide a complete set 
of coordinates to describe this subspace. We will use the 
theorems of Karhunen-Loeve and Mercer and will limit the 
demonstration to the dc waveform for simplicity and without loss 
of generality. We start by considering the expansion of a random 

20 complex sequence {z(n), Vn} in a series of waveform coordinates 
consisting of time translates of the mother waveform y/ . The 

sequence {z(n), Vn} is a zero-mean stationary random process 
which is orthonormal over the sample interval n M" and has a 
frequency spectrum which is flat and extends over the frequency 
25 range 1/LT which is centered at baseband corresponding to a zero 
frequency. This means the {z(n), Vn} cover the subspace of V 
corresponding to the scale of our new dc waveform y/ and its time 

translates { y/(n -qM) = y/ q {n\ Vq } . In addition, the V is now 
considered to be extended over a time interval which is 
30 relatively large compared to the N-dimensional t-f space in FIG. 
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1,2 to avoid end-effects on the analysis. We start by 
approximating the sequence {z(n), Vn) by the (z(«),Vn) where: 



z(n)= £Z q ^(n-qM) 
q 



10 



= £Z q + (n,-q)M) 

q 

where the complex coefficien ts {z q }are 

derived from the original sequence using the ( 12 ) 

orthonormality properties of y/ 
Z q = Zz(n) ^*(n-qM) 

n 

= Zz(n)^ q (n) 

n 

The following equations prove that the coefficients {Z q , Vq} are 
orthonormal : 



Z q Z q ,* = Sz(An-qM)z*(An-q f M)^ q (Any q , *(An) (13) 

An 

= 8<tf Z V q ( n )| 2 since the sequence {z} is orthonormal 

0 

for qM time translates 
= £ qq . with normalization of the energy of y/ 



Equations (12) and (13) together prove the Karhunen- 
15 Loeve's theorem which proves the following equation for the 
accuracy in approximating the stochastic sequence {z(n), Vn} by 
{z(n),Vn}. This accuracy is expressed by the expected "E(o)" 
squared error " (o) " in this approximation: 



20 
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E[z(n) - z(«)| 2 } = 1 - J>(«o + n x Mf 



(14) 



We need to prove that the right hand side of this equation is 
zero which then proves that the approximating sequence is equal 
to the original sequence in the mean-square sense. In turn this 
proves that the new waveform coordinates { y/ q , Vq } are a basis for 
the original sequence {z(n), Vn} which is our goal. 

The right hand side of equation (14) when set equal to zero 
expresses Mercer's theorem so our goal is to prove Mercer's 
theorem. To do this we use the DFT of y/ in equation (11) and the 
coordinates in (10) to evaluate the right hand side of equation 
(14) . We find 




where T = <VLf£Wi 



= sin(;r Ak 0 )/Isin(;r Ak 0 1 L) 
= 1 for Ak 0 = 0 



(15) 



= 0 otherwise 
This proves that 




o 



which proves that equation (1S14) reduces to 



E[z(n)-z(n)\ 2 \ 



= 0 



(16) 
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which as per the above proves that the set of multi-resolution 
waveforms is a basis. This proof easily generalizes to the 
multi-resolution waveforms at all of the scales {p} and time 
translates {q}, and to the expansion of the harmonic design 
5 coordinates over 2L, 3L, ... as required by the application. 



New — Our new waveforms have multi-scale properties that we 
will demonstrate. We start with the observation that the Fourier 
domain design for the new waveforms provides a natural and easy 

10 way to derive the complete set of waveforms { y pq , Vp,q } for the 
space V from the design of the dc waveform y/ by using the same 
invariant set of Fourier domain harmonic design coordinates 
{ x ¥ k ,Vk 0 } derived for the dc multi-resolution waveform. This 
demonstration requires that we show 1) how the multi-scale 

15 transformations are implemented with the design in the Fourier 
domain, and 2) how the waveform design remains invariant under 
scale changes . 

First consider the multi-scale transformation which derives 
20 the waveforms at the scale and shift parameters M p,q" for from 
the dc waveform at scale M p=0" and centered at the origin "q=0". 
We begin by extending the parameters and coordinates in equation 
(10) , to include both scaling and subsampling or decimation in a 
form that is equivalent to the iterated filter bank construction 
25 which is used to derive ^t^e — current Wavelet waveforms using the 
filter scaling functions. Starting with the coordinates at scale 
«p=0" the parameters and coordinates at scale "p=p" are given by 
the equations: 

30 



.pp=0 parameters and coordinates 

n = n 0 + ni M (17) 

n 0 = a 0 + ai 2 + ... + a m -i 2 m_1 
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= M points 
M = 2 m 

ni = b 0 + bi 2 + ... +bi-i 2 1 " 1 

= L points spaced at M sample intervals 



p=p parameters and coordinates 

2' p n(l2 p ) = scaled by "2" p " 

and subsampled or decimated by 2 P :1 
10 = n 0 (p) + ni(p)M 

2" p n 0 {l2 p ) = scaled by "2" p " 

and subsampled or decimated by 2 P :1 
= n 0 (p) 

= a p + ap~+i2 + ... + a p+m -i 2 P 
15 = M points spaced at 2 P sample intervals 

2" p n x (l2 p ) = scaled by "2" p " 

and subsampled or decimated by 2 P :1 
= ni(p) 

= bp + b p+1 2 + ... + b p+1 -! 2 P+1_1 
20 = L points spaced at M2 P sample intervals 



together with the observation that the sampling interval "T" is 
increased to "2 P T" under the scale change from "p=0" to "p=p" and 
25 subsampling or decimation from "1:1" to "2 P :1". Combining these 
equations with the analytical formulation in (7) and the Fourier 
domain representation in (11) enables the waveforms at the 
parameters "p,q" to be written as a function of the Fourier 
domain harmonic design coordinates: 

30 
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= 2" ,V2 K2" p n-qM)e i2rfc(p ' r)nT 
= (2~ p/2 /N') W£ (aWH|M> e 



i2*f c (p,r)n(p)2 p T 



(18) 



for all admissible scale, translation, and frequency index 
parameters "p, q, r" . 

Next we need to demonstrate that the frequency domain 
design in (11) remains invariant for all parameter changes and in 
particular for all scale changes. This multi-scale property 
expresses the accordion behavior of the design in that the 
Wavelets at different scales are simply the stretched and 
compressed versions of the mother waveform with the appropriate 
frequency translation indices. This multi-scale invariancy means 
that the design for a M=16 channel filter bank remains the same 
for M=100 or 14=10,000 channel filter banks, when the overlap L 
and the performance goals remain constant. To demonstrate this 
invariant property across scales, we consider the mag— multi- 
resolution waveform at scale "p" with the other parameters set 
equal to zero for convenience "q=0, r=0" and without loss of 
generality. The Fourier domain frequency response *P(/) can be 
evaluated starting with the original formulation in equation (7) : 



DFT 



at "p, q=0, r=0 



(19) 
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k 0 n(p) 



-UN*2>T -k 0 )n(p) 



S inQ (fN'2 p T-k 0 )) 
N'sin(;r (fN'2 p T-k 0 )/N f ) 



= 2] ^ [Harmonic interpolation for " k 0 ft ] 



k 0 



This only differs from the harmonic representation in equation 
5 (9) in the restriction of the design coordinates to the subset of 
harmonic coefficients ,Vk 0 } and the stretching of the time 

interval to "2 P T" corresponding to the scale "p=p". The harmonic 
interpolation functions are observed to remain invariant over 
scale changes upon observing that the frequency scales as 

10 "f~l/2 p T'' which means the "f requency*time" product remains 
invariant with scale changes as per the fundamental property of 
the waveforms. This means the harmonic interpolation functions 
remain invariant with scale change and therefore the frequency 
response remains an invariant. This demonstrates the waveform 

15 design is an invariant across the waveform scales which means we 
only need a single design for all scales or resolutions of 
interest. 

LS design algorithms for new waveform will be described to 
20 illustrate the advantages our ftew — waveform has over current 
designs. The two LS algorithms described are the eigenvalue and 
the gradient search which respectively can be reduced to 

algorithms which are equivalent to the original current 

eigenvalue [7] and Remez-exchange [8] waveform design algorithms 
25 for application to a uniform filter bank. We consider the t-f 
space which is spanned by a uniform polyphase filter bank 
consisting of M channels at the frequency spacing f w =l/MT where T 
is the digital sampling interval, and the filter waveform FIR 

24 



time response is stretched over L sampling time intervals T s . 
This polyphase filter bank is ideally decimated which means the 
filter ouput sample rate 1/T S is equal to the channel-to-channel 
spacing 1/T 3 =1/MT, equivalent to stating that there is no excess 
5 bandwidth a=0 . Our design for this topology is immediately 
applicable to an arbitrary set of multi-resolution filters 
through the scaling equation (18) which gives the design of our 
waveform at arbitrary scales in terms of our design of the dc 
waveform. 

10 

For this polyphase filter bank used to construct the dc 
waveform or filter impulse response, our LS example design 
algorithms will use 5 metrics consisting of the 2 prior art 
passband and stopand metrics, and the 3 new metrics consisting 

15 of the ISI, ACI, and QMF, and solve the LS minimization problem 
using as design coordinates the subset of harmonic coordinates 
which are a basis. Since our 2 — two example LS design algorithms 
only differ in the use of an eigenvalue LS optimization and the 
use of a gradient search LS optimization, the flow diagrams for 

20 the construction of the cost functions and the solution for the 
optimal waveform will be identical. However, there are 
differences in the construction of the cost functions from the 
respective metrics and in the iterative solution mathematics. 
Both LS solutions for the harmonic design coordinates minimize 

25 the weighted sum of the error residuals or cost functions from 
the 5 metrics. These design coordinates — are the Fourier 
harmonics {^,Vk 0 } for {k 0 = 0, 1 , L-l } - Resulting algorithms 
are easily extended to the applications requiring the design 
coordinates to cover 2L, 3L, ... harmonics. 

30 

Frequency domain design coordinates are related to the 
waveform time domain digital samples or coordinates as follows. 

Mappings of timo — £3— frequency to time (20) 

25 



Time domain design coordinates {y/(n),\/n} are real and 

symmetric and can be represented by the reduced set {h t (n), 
n=0,l,...,ML/2} 

h t (n) = y(0) for n=0 

= 2y/(n) for n=l, 2 / ... / ML/2 

= time domain design coordinates 

Frequency domain harmonic design coordinates {^ ? Vk 0 } are 

real and symmetric and can be represented by the reduced 
set {h f (k) , k=0,l,,...,L-l} 

h f (k) = y/ ko for k=k o =0 

= 2y/ kQ for k=k 0 =l,2, ...,L-1 

= frequency domain design coordinates 

Mapping of the frequency coordinates {h f (k), k=0 , 1 , L-l } 
into the time coordinates {h t (n), n=0, 1, ... / ML/2} is defined 
by the matrix transformation to within a scale factor 
h t = B h f 

where 

h f = (h f (0), ...,h f (ML/2) ) t transpose of column vector 
h t = (h t (0), ...,h t (ML/2) ) t transpose of column vector 
B = (ML/2 +1) x L matrix 

= [ B kn+ i /k+ i ] matrix of row k and column n 

elements B kn 

B* n+ i,k + i = 1 /ML for n=l 

= 2 cos(27ikn/ ML) otherwise 



Mappings — of the frequency coordinates — R* *(k) / — k-0, 1, L - l } 
into — the — time — coordinatco — H Mn) / — n-0, l,...,ML/2] — aftd — vice 
versa aro defined by the matrix transformations 
h*__ =-B — h# 
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where 

(h*(0), ...,h±(ML/2) ) fe — tranopooc of column vector 

.=—+&t ( 0 ) , ...,h fe (ML/2) )* — tranopooc of column vector 

B — = — (ML/2 — I 1 ) x L matrix 

= [ B ^H matrix of row k and column n clcmcnto B^ 



-B ^ 1 /ML + 1 for n-1 

- 2 coo (2rc]cn/ ML + 1) othcrwioc 



B" 1 ~ L x (ML/2 + 1) matrix 

10 = — (rB-^l matrix of row k and column n clcmcnto - B^j- 

for n-1 



2 coo (2rckn/ ML f 1 ) othcrwioc 



and 



B~* ~B - I L x L identify matrix 

15 B— B" 1 ~ I (ML/2 i l) x (ML/2 i l) identity matrix 



wherein N' =ML+1 has been replaced by ML since a single end point 
has been added to the FIR to make it symmetrical for ease of 
implementation . 
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Passband and stopband metrics and cost functions are 
derived with the aid of FIG. 3 which defines the power opcctral 
denoity (-PSDf parameters of interest for the passband and 

25 stopband of the PSD ^(a?) for communications applications. 
Requirements for radar applications include these listed for 
communications. Referring to FIG. 3 the passband 11 of the 
waveform PSD is centered at dc (f=0) since we are designing the 
dc or baseband waveform, and extends over the frequency range ©p 

30 extending from -co p /2 to +©p/2 12 in units of the radian 

frequency variable o=27CfT 13 where T io the digital sampling 

interval defined in FIG. 1 . The frequency space extends over the 
range of f=-l/2T to f=+l/2T which is the frequency range in 
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FIG 1 translated by -1/2T so that the dc waveform is at the 
center of the frequency band. Quality of the PSD over the 
passband is expressed by the passband ripple 14. Stopband 15 
starts at the edge 16 of the passbands of the adjacent channels 
5 + /_ 0)a /2 16 and extends to the edge of the frequency band ©=+/- 
7i 17 respectively. Stopband attenuation 18 at +/-co a /2 
measures the PSD isolation between the edge of the passband for 
the dc waveform and the start of the passband for the adjacent 
channels centered at + /-g> s 19- Rolloff 20 of the stopband is 

10 required to mitigate the spillover of the channels other than the 
adjacent channels, onto the dc channel. Deadband or transition 
band 21 is the interval between the passbands of contiguous 
channels, and is illustrated in FIG. 3 by the interval from (Dp/2 
to (o a /2 between the dc channel and adjacent channel at © a . 

15 Waveform sample rate <o s 22 is the waveform repetition rate. 
For the LS example algorithms, the waveform sample rate is equal 
to the channel-to-channel spacing for zero excess bandwidth. 
Therefore, 1/T S = © s /27iT=l/MT which can be solved to give co s =27i/M 
for the radian frequency sampling rate of the filter bank which 

20 is identical to the waveform repetition rate. 

We start by rewriting the DFT equations for the dc waveform 
in (11) as a function of the {h f (k), k=0, 1,...,L-1} 

y/(o>) = X^(n)cos(rutf) 

n 

25 = c' B h f using (20) and the definition of the vector "c" (21 ) 

c = ( 1 , cos(<y), cos( (ML/2) co) )' transpose of column vector 

which is equivalent to the equation for if/((o) in (9) expressed as 
a linear function of the {h f (k), k=0, 1,..., L-l} . However, this 
30 functional form suggested by the eigenvalue formulation [7] is 
more convenient to analyze. An ideal w c" vector "c r " will be 
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introduced for the passband and the stopband in FIG. 3-r in order 
to identify the error residual Sy/(6)) at the frequency "G)" in 
meeting the ideal passband and stopband requirements. The ideal 
PSD is flat and equal to "1" for the passband, and equal to "0" 
for the stopband. We find 



Error residuals for passband and stopband (22) 

c r = (1,1,..., l) fc passband ideal "c" 
{0,0,.. .,0)* stopband ideal "c" 
Sc = c r - c error vector 

8y/(<D) = Sc l Bh f 

= Residual error in meeting the ideal spectrum at n a>" 



The LS metric for the passband and stopband can now be 
constructed as follows for the eigenvalue and the LS optimization 
(or cquivalcntly, the LS algorithm) design algorithms- 



Passband and stopband metrics (23) 



J (band) = — - — J|^(a;)| dco Eigenvalue 



band ^ 



= h' f R h f Eigenvalue 
= ||<V ||- LS 



—where 

band = [0, co p ) passband 
= (co s , 7i] stopband 



(23) — continued 
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1 



R 



band 



bond 



L x L matrix 



8i(/ 



( 8\|/(g>i) 8vj/ (o) u ) ) t 

vector of error residuals at the 



frequencies o>i , 



... , 



G) u across the band 



IN 



norm or length of the vector (o) and which 



includes a cost function for the errors of 



the individual components 



where it is observed that the eigenvalue approach requires that 
the LS metrics be given as quadratic forms in the design 
coordinates {h f (k), k=0 , 1 , L-l } whereas with the LS approach it 
is sufficient to give the LS metrics as squared vector norms 
with imbedded cost functions or an equivalent formulation. 

QMF metrics express the requirements on the deadband that 
the PSD' s from the contiguous channels in FIG. 3 add to unity 
across the deadband [©p , o s ] in order that the filters be QMF 
filters. By suitable modification of the error vector 8c, the 
previous construction of the passband and stopband metrics can be 
modified to apply to the deadband. This is a Wavelet requirement 
on the frequency coverage of the Wavelet basis— as — observed in 
Table 1. . We find 



Deadband metrics (24) 



J (deadband) 




deadband 



Eigenvalue 



deadband 



30 



(24) continued 
= h f Eigenvalue 

= ||<V ||- - LS 

where 

8c = c r - c (co) - c(7i/M - co) 



where c(o)=c as defined in (22) and (23), and c(ti/M - ©) = c at 
the offset frequency "7t/M - <o" corresponding to the overlap of 
10 the contiguous filters over the deadband. 

Orthonormality metrics measure how close we are able to 
designing the set of waveforms to be orthonormal over the t-f 
space, with the closeness given by the ISI and the ACI. ISI is 

15 the — non - orthogonality — error — between — channel — output — samples 

separated by multiples of the — sampling interval 1/MT seconds 

where T is the sample time and M is the interval of contiguous 
samples - — AC^ — i-s — the — non - orthogonality — error — between — between 
channel — output — samples — within — a — channel — aad — the — samples — if* 

20 adjacent — channels at the — same sample time and at sample times 

separated by multiples of the sample interval. As observed as 

noise contributions within each sample in a given channel, the 

ISI is the noise contribution due to the other received waveforms 
at the different timing offsets corresponding to multiples of the 

25 sampling interval, Likewise, the ACI is the noise contribution 

due — fee — the — other — waveforms — ±rt — adjacent — channels — at — the — same 
sampling time and at multiples of the sampling interval. 
ISI and ACI errors are fundamentally caused by different 
mechanisms and therefore have separate metrics and weights to 

30 specify their relative importance to the overall sum of the LS 
metrics. ISI is a measure of the non-orthogonality between the 
stream of waveforms within a channel as per the construction in 
FIG. 3. On the other hand, ACI is a measure of the non- 
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orthogonality between the waveform within a channel and the other 
waveforms in adjacent channels. This means the stopband 
performance metric has a significant impact on the ACI due to the 

sharp rolloff in frequency of the adjacent channel 7 and the 

5 ACI metric io then a measure of the residual — non - orthogonality 
due to the inability of the otopband rolloff in frequency from 
completely eliminating the ACI errors. 

We assume that the received waveform is identical to the 
10 filter waveforms and is transmitted at the filter output sample 
intervals equal to ±/MT seconds. The second assumption means we 
are assuming the receiver is synchronized with the received 
signal. Since — there — — J*e — information — loot — by — sampling 
aoynchronouoly with the — received wavcf orm, — we — enee — free to make 
15 this synchronization assumption without loss of generality. — ISI 
metrics are derived in the following set of equations. 

ISI metrics (25) 

20 Mapping of h f into i|/ 

x ¥' = (i|/(-ML/2), ... , \|/(ML/2)) t transpose of 

column vector ¥ 

= H h f 

H = (ML+1) x L matrix of elements ft^H mfk+ i 

25 H m , k+ i = 1 for ftk=4-0 

= 0.5 cos (27tkn/ML+l) otherwise 
where m = n + ML/ 2 + 1 
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Offset matrix A 
A = L x (ML+1) 
A k+1 ,™ = [0 0 ... 0 



matrix of elements A k+ i /ft - 
0 ] row vector k=0-9- 
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= [0 —0 — ... — y (-ML/2)—. V|/ ( (L-k)M+l) ] row vector k+1 

(25) — continued 

12 kM LM+1 column m element 

of row vector 

5 

ISI error vector 5E 
8E = L x 1 column vector 
= AH h f 

10 ISI metric 

J (ISI) = 8E fc 5E Eigenvalue 

= Non-linear quadratic function of h f 

= WseU ls 
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ACI metrics are derived using the ISI metric equations with the 
following modifications. 

ACI metrics 
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Mapping of h f into \\) is the same as developed for ISI 
Offset matrix A elements are changed as follows to apply to 
Channel 1: 

= [0 0 0 0 ] — row vector k=4-0 

25 =[0 0 ... h/(-ML/2)^ ... i[/((L-k)M+l) W£~ k) ] row 

vector k+1 

1 2 kM+1 LM+1 column m element 
of row vector 

which means the ACI error vector 8E is 
30 SE = L x 1 column vector 

= A H h f 

ACI metric for the two contiguous channels 
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J(ISI) = 2 SE 1 8E Eigenvalue 

= Non-linear quadratic function of h f 



= 2 \\SE\\ 2 - LS 



5 where the factor "2" takes into account there are two contiguous 
channels or one on either side of the reference channel 0 in FIG. 
3. Because of the fast rolloff of the frequency spectrum the 
addition of more channels into the ACI metric is not considered 
necessary, although the functional form of the ACI metric in (26) 
10 allows an obvious extension to any number of adjacent channels 
which could contribute to the ACI. 

Cost function J for the LS algorithms is the weighted sum 
of the LS metrics derived in (23) , (24) , (25) , (26) . The LS 
15 algorithms minimize J by selecting the optimal set of frequency 
coordinates {h f (k),Vk} for the selected set of parameters used 
to specify the characteristics of the dc waveform, frequency 
design coordinates, LS metrics, and weights. Cost function and 
optimization techniques are given by the equations 
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Cost function J 



(27) 



J 




= weighted sum of the LS metrics J (metrics) 
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where 



metrics = passband, stopband, deadband, ISI, ACI 
{w (metrics) , ^metrics} = set of weights 



^{metrics) - 1 normalization 



X 
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Optimization goal 

Goal: minimize J with respect to the selection of 



the {h f (k),Vlc} 
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Optimization algorithms 

Two algorithms are the Eigenvalue and the LS optimization 
5 where the eigenvalue optimization algorithm uses the non- 

linear quadratic formulations of the LS metrics and the LS 
optimization algorithm uses the squared norm formulations 
for the LS metrics. 

10 

FIG. 4 is a summary of the LS metrics and the construction 
of the cost function J. Design parameters 23 arc the input and 

output design paramctoro. Input parameters are L, the number of 

polyphase — channels — M,_ — cquivalcntly — fcke — number — e-f — digital 

15 samples at spacing T over the oymbol interval T e =M^4 the length 

of the FIR time response for the waveform in units of L which are 
the number of digital samples per waveform repetition interval T e 
so that the total number of digital samples for the symmetric FIR 
time response is equal to N'=M L+1, number of DFT samples per FIR 

20 length n_fft for implementation of the LS algorithms, passband 
radian frequency stopband radiam frequency CQ a , waveform 

repetition rate in radian frequency oo s , selection of the set of 
design coordinates {h f } to be used in the optimization, and the 
metric weights {w (metrics) } . Output parameters are the set of 

25 harmonic design coordinates {h f } that minimize J. Band metrics 

24 are the passband, stopband, and deadband metrics defined in 
equations (23), (23), (24) respectively. Interference metrics 

25 are the ISI and ACI metrics defined in equations (25) and 
(26) respectively. LS cost function J 26 is the weighted 

30 linear sum of the metrics defined for the band 24 and the 
interference 25 as defined in equation (27) . 

FIG. 5 is a flow diagram of the LS recursive solution 
algorithm. There are two loops with the topology constructed so 
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that the outer loop 27 is an iteration over the set of metric 
weights {w (metrics )} , and the inner nested loop 28 is the 
recursive or iteratived LS solution to find the optimal {h f } for 
the given design parameters and weights where each step refers to 
5 an iteration step in the solution. A recursive LS solution is 
required to be able to solve the highly non-linear interference 
metrics as well as the band metrics when one chooses to use non- 
linear techniques to construct these metrics. 

10 



15 

Consider the inner nested loop 28. The recursive LS 
solution starts with the initial step i=0 29 which begins with 
the selection of the input design parameters and weights and the 

20 selection of an initial set of values for the {h f } 30. Next the 
band metrics are calculated 31. Since this is the initial step 
i=0 32 the cost function J is restricted to the linear band 
metrics in order to fine the linear approximation to {h f } 33 to 
initialize the non-linear solution starting with step i=l 34 

25 wherein the highly non-linear interference metrics are included 
in J. Following the signal flow, for step i=l the band metrics 
31 and interference metrics 35 are calculated, weighted, and 
summed to form the cost function J 36. An iterative solution 
algorithm 37 finds the best approximation in step i=l to the 

30 {hf} which minimizes J using the approximation to {h f } from the 
previous step i=0 to linearize the search algorithm coefficients 
for the eigenvalue and LS step i=l iteration. If there is no 
convergence to the correct solution for {h f } this recursive 
sequence of calculations is repeated for step i=2 34. This 

35 recursive solution technique is repeated for subsequent steps 



until there is convergence 38 whereupon one exits this inner 
loop . 

Consider the outer loop. After exiting the inner loop 38 
5 the solution for {h f } is tested to see if it meets the 
performance goals 39. If not, a new set of metric weights is 
selected and the inner loop is initialized i=0 40 and the 
inner loop is used to find the next solution for {h f } . This 
process continues until the performance goals are met or 
10 adequated approximated 39 whereupon the algorithm is exited 
with the final solution set {h f } 41. 

Applications of this new invention to both communications 
and radar will be given using these example algorithms and other 

15 algorithms supported by this invention. These new waveforms are 
considered for the applications: l)to replace the square-root 
raised cosine waveform (sq-root rc) which is extensively used for 
the third generation (3G) CDMA communications, 2) to replace the 
Gaussian minimum shift keying (GMSK) waveform for constant 

20 amplitude bandwidth efficient (BEM) applications, and 3) as a 
candidate waveform for synthetic aperture radar (SAR) and real 
aperture radar (RAR) applications. 

CDMA communications applications for the current and the 
25 new 3G CDMA considers the use of a waveform dcoigncd with this 
ftew — invention — as — a — possible — replacement — — the — sq-rt rc 
waveforms with bandwidth expansion parameter a = 0.22 to a - 0.40 . 
This notation means that for a = 0.22 the spectral efficiency is 
(symbol rate/bandwidth) = 1/1 + a =1/1.22 = 0.82 = 82%. A basic 
30 advantage of t&e -our waveform for 3G CDMA is the potential for a 
symbol rate increase within the same bandwidth with an increase 
in the spectral efficiency to «100% depending on the application 
and operational constraints. The dc power spectral density or 
power spectrum (PSD) of the waveform is compared to the PSD for 
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the sq-rt r-c in FIG. 6. Plotted are the measured PSD in dB 42 
versus the frequency offset from dc expressed in units of the 
symbol rate 43 Plotted against the normalized frequency offset 
are the dc PSD for the new waveform 44, the sq-rt r-c with 
5 a=0.22 45, and the sq-rt r-c with a=0.40 46. It is observed 
that the PSD for the new waveform rolls off faster that that for 
the sq-rt r-c which means that the new hh? — waveform will 
oupport of fers the potential for an increased symbol rate for a 
given available frequency band^ while — oatiof ying — fc£e — inherent 
10 rcquircmcnto for low ISI and MAI — (multiple acceoo intcrf Grcncc) . 

Constant amplitude BEM application of the new waveform 
indicates that it is a viable candidate for replacing the current 
preferred modulation waveform which is the GMSK. The GMSK finds 

15 applications for transmitters which operate their HPA(s) 
amplifiers in a saturation mode in order to maximize their 
radiated power from the HPA(s), and which require a BEM PSD to 
avoid excessive spreading of the transmitted power. Simulation 
data for the dc PSD is plotted in FIG. 7 for the new waveform BEM 

20 and the GMSK. Plotted are the measured dc PSD in dB 47 versus 
the frequency offset from dc expressed in units of the bit rate 

48- Plotted against the normalized frequency offset arc the dc 

P&B for the new waveform BEM 49 and the GMSK 50- for a length 
parameter L=10 where L is the length of the phase pulse in terms 

25 of the phase pulse repetition rate. Tfee — oignif icancc — — this 
This example data — is — that — &he — aew — waveform — has — indicates the 
potential to bo — designed to of our waveform to offer a PSD which 
is less spread out than the current GMSK, and therefore an 
improved BEM waveform. 

30 

Radar RAR and SAR application of the waveform indicates 
that it is a viable candidate to replace the current chirp 
waveforms for wideband signal transmission, when combined with 
pseudo-random phase codes. Results of the simulation for the new 
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waveform and an unweighted frequency chirp waveform are given in 
FIG. 8. Plotted are the ambiguity function for the new waveform 
51 and the unweighted frequency chirp waveform 52. The dc 2- 
dimensional radar ambiguity function 53 is plotted as a function 
5 of the frequency offset in units of fT p and the time offset in 
units of t/T c where T p is the phase-coded radar pulse length or 
length of the phase code and T c is the phase code chip length. 
The chip length is identical to the waveform repetition interval 
T s so that T c = T s . it — is — obocrvcd — that — th eThis example 
10 indicates our new waveform has the potential for significant 
improvements in the ambiguity function and by implication in the 
performance . 

Preferred embodiments in the previous description is — are 
15 provided to enable any person skilled in the art to make or use 
the present invention. The various modifications to these 
embodiments will be readily apparent to those skilled in the art, 
and the generic principles defined herein may be applied to other 
embodiments without the use of the inventive faculty. Thus, the 
20 present invention is not intended to be limited to the 
embodiments shown herein tetand— is not to be accorded the wider 
scope consistent with the principles and novel features disclosed 
herein. 

25 
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DRAWINGS AND PERFORMANCE DATA 



PIG. 1 Wavelet Tiling of an N- Point Digital t-f Space 
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FIG. 2 Wavelet Iterated Filter Bank for Tiling t-f 
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FIG. 3 PSD Requirements for Communications 
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FIG, 4 



LS Metrics and Cost Function 
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Design Parameters 




Passband Eq. 23 
Stopband Eq. 23 
Deadband Eq. 24 



ISI Eq. 25 
ACI Eq. 26 



LS Cost equation 
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FIG. 5 LS Recursive Solution Algorithm 
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FIG. 6 PSD fox New Waveform and Square-Root Raised-Cosine 
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FIG. 7 PSD for New Waveform BEM and GMSR 




FIG. 8 Radar Ambiguity Functions of New Waveform 
and Unweighted Chirp Waveform 



New Wavelet Ambiguity Function 51 Unweighted Chirp Ambiguity Function 52 
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